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Abstract 



Recent applications in queuing theory and statistical mechanics have isolated the process formed 
by the eigenvalues of successive minors of the GUE. Analogous eigenvalue processes, formed in 
general from the eigenvalues of nested sequences of matrices resulting from random corank 1 
projections of classical random matrix ensembles, are identified for the LUE and JUE. The 
correlations for all these processes can be computed in a unified way. The resulting expressions 
can then be analyzed in various scaling limits. At the soft edge, with the rank of the minors 
differing by an amount proportional to N'^^^, the scaled correlations coincide with those known 
from the soft edge scaling of the Dyson Brownian motion model. 



1 Introduction 



Since the pioneering days of Wigner, Gaudin, Mehta and Dyson (see Porter \31\ for a collection of 
papers from this period), random matrix theory has shown itself to be perhaps the richest source 
of exact solutions for correlation and distribution functions of all statistical mechanical models. 
Many of the discoveries of this type have their motivation in new applications of random matrix 
theory, unknown in the pioneering days. A case in point is the recent work of Johansson and 
Nordenstam |22[ [27] , who compute the exact form of the correlation functions for the coupled 
eigenvalue sequences obtained from the principal minors of Gaussian unitary ensemble (GUE) 
matrices. 

The motivation for studying this GUE minor process begins with a work of Baryshnikov [1]. 
Some years earlier Glynn and Whitt [15j had studied the problem of computing the distribution 
of exit times from a queueing system, in the limit the number of queues tends to infinity but 
the number of jobs remains finite. It was proved that for general i.i.d. service times the scaled 
distributions of the exit time of the kth customer from the final queue could be written as 



Here each Bi denotes an independent standard Brownian motion, while the condition X is that 



By studying the particular case of exponential waiting times Baryshnikov was able to show 
that {-Dfc}/c=i,2,... could alternatively be specified as the joint distribution of {//fc}, where //fc is 
the largest eigenvalue of the kth principal minor of an infinite GUE matrix X with probability 
density function (PDF) proportional to exp(— X^/2). Johansson and Nordenstam [22] give other 
occurrences in statistical mechanics of essentially the same process relating to the eigenvalues of 
minors of GUE matrices (referred to as the GUE minor process). These are in the specification 
of certain point processes relating to the boundary region (neighbourhood of the frozen zone) of 
random domino tilings of the Aztec diamond [20], and to random lozenge-tilings of a hexagon 
[21] . In the equivalent language of stepped surfaces, Okounkov and Reshetikhin [29] make similar 
observations. The recent work of Borodin, Ferrari and Sasamoto [5] encounters this process in 
the context of studying the dynamics of the asymmetric exclusion process. 

The broader setting of the measures encountered in the studies of the above statistical 
mechanical models relates to the Robinson-Schensted-Knuth (RSK) correspondence (see Section 
12. ip below. Consideration of the structures inherent therein [21 |T4] identifies natural extensions 
of the GUE minor process. One such class of extensions replaces the Gaussian weight in the 
latter by the Laguerre or Jacobi weights, which can all be realized by a sequence of projections 
onto random complex hyperplanes, giving rise to so called classical projection processes (the 
Gaussian, Laguerre and Jacobi weights are all classical from the viewpoint of the theory of 
orthogonal polynomials). Our main point in the present paper is that the correlations for the 
classical projection processes can be computed exactly in a unified way. The essential ingredient 
here is the Rodrigues formula for classical orthogonal polynomials. Moreover, known asymptotic 
formulas for the latter allow for the evaluation of scaling limits of the correlations. 

We begin in Section 2 by recalling how the RSK correspondence relates to a certain statistical 
mechanical model of last passage times. In Section 3 the occurrence of special cases of the joint 
PDF in some random matrix setting is noted. The correlations for the classical projection 
process are calculated in Section 4, and their scaling limits analyzed in Section 5. 
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2 A joint probability density associated with RSK 



2.1 The case of general parameters 

The Robinson-Schensted-Knuth (RSK) correspondence gives a bijection between ni x n2 non- 
negative integer matrices [xij] (rows countered from the bottom) with entry (ij) weighted 
{aibjY^'^ , and pairs of weighted semi-standard tableaux (weights {aj}, {hj}) of shape n = 
(//I, . . . ,//n))- Results from [ISl [2] associate a probabilistic model to the RSK correspondence. 
Identify with each lattice site {ij) a random non-negative integer variable xij chosen from the 
geometric distribution with parameter aibj, so that 

Pr(x,j =k) = {l- a,hj){a,hjf. (2.1) 

For given non-negative parameters {oj}, and given ni,n2 G Z+, a probabilistic quantity of 
interest is the sequence of last passage times 

L^^\ni,n2) = max Xij, / = 1, . . . , max(ni, n2). (2-2) 

(rd*)' 

Here (rd*)^ denotes the set of / disjoint (no common lattice points) rd* lattice paths, which 
in turn are defined as either a single point, or points connected by segments formed out of 
arbitrary positive integer multiples of steps to the right and steps up in the rectangle 1 < i < ni, 
1 < J < n2. 

A crucial feature of the RSK correspondence is that the length of the first row of the semi- 
standard tableaux pair is equal to (|2.2p in the case I = 1, and thus /ii = L^^\ni,n2)- More 
generally, all row lengths are determined by (j2.2p according to [TT] 

W = L«(ni,n2)-L('-i)(ni,n2) (2.3) 

with L('')(ni,n2) := 0. Thus the distribution of the last passage times is fully determined by 
the distribution of {/i;}. 

A well known fact relating to the RSK correspondence is that the probability an ni x n2 non- 
negative matrix with elements chosen according to (j2.ip corresponds to a pair of semi-standard 
tableaux with shape /i, one of content ni, the other of content n2, is given by |24j 

ni n2 

P{^J) = n ~ aihj)s^{ai, a„Js^(6i, . . . , hn^), (2.4) 
i=ij=i 

where denotes the Schur polynomial. A lesser known fact is that for ni > n2 the joint 
probability that an ni x (n2 + 1) non- negative integer matrix with elements chosen according 
to (j2.ip corresponds to a pair of semi-standard tableaux with shape ^, content ni and n2 + 1, 
and the bottom left sub-matrix corresponds to a pair of semi-standard tableaux with shape k, 
content ni and n2 is |14j 

Y\W_{'^-aihj)s^,{ai,...,ani)s,,(hi,...,hn^)h^^'^l ' ' > n) (2.5) 

i=i j=i 

where, with x(^) = 1 if A is true, x(^) = otherwise, 

x(/i > k) := xifJ-i > > /"2 > > • • • > /^m > > 0). (2.6) 
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It follows from (j2.4p and ()2.5p that given the pair of semi-standard tableaux corresponding to 
an ni x 722 matrix (ni > 71-2) has shape k, the probability of the pair of semi-standard tableaux 
corresponding to the ni x (77-2 + 1) matrix, obtained by adding an extra row to the existing 
matrix, having shape fi is 

\ / ^ \TT/-i u \ . . . , a„J E"ii(Mj-'«j)+/^n2+i v\ 

''k V*^! ) • • • ) '^ni ) 

Let us now seek the joint probability that with ni > n2 + p, an ni x (712 + ^j) non-negative 
integer matrix with elements chosen according to (j2.ip is such that the principal rii x (n2 + s) 
sub-blocks (s = 0, 1, . . . ,p) correspond to pairs of semi-standard tableaux with shape fi^^\ This 
is computed from (j2.4p and (|2.7p according to 

s=0 



ni n2+P 

i=l j=l 

P s-^n2+3-l, (s) I (s) 

xnfe ^ '"^■^^x(m(^)>/x(^-^)). (2.8) 

s=l 



2.2 Specializing the parameters 

In [13] the joint probability (|2.5p is specialized to the case of a geometrical progression of 
parameters 

(ai,...,a„J = zt^, . . . 

(61,..., = (2.9) 

This is a preliminary step for taking the so called Jacobi limit, in which a joint probability of a 
type known from random matrix theory is obtained. In this subsection the parameters will be 
specialized according to ()2.9p for the more general joint probability (|2.8p . 

Now for (j2.8p to be non-zero we require £(/i^P^) < n2+p- Under this circumstance, we deduce 
from [14, Eq. (2.26)] with k ^ fJ-^^K ^2 1— > ni, ni ^ n2 + p, rj 1— > h^^ := /u^-^"* + n2 + p — j that 

+-E"^0-l)(n2+p-i) 

S„(rt(l,t,...,t"^-l) = r^^=i J0-V(n2+P)EA ' . 

11/=1 (*>*)« 

i=l j=l ^ ' '^/i-*'' i<i 



This makes explicit the first Schur polynomial factor in (|2.8p . In regards to the second Schur 



polynomial factor, making use of [141 Eq. (2.18)] with n 722, n* 1— > 722, A 1-^ fJ'^^K hj h^^^ 



(0) , 

l^j +n2- 3 gives 



.-E;ii(i-i)("2-i) "2 



Use of these results in (12. 8p allows the following results to be deduced. 
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Proposition 1. Let ni > n2 + p- On each site of the ni x (^2 + p) square lattice specify a 
random non-negative integer Xij according to the specification 

Pj:{xij = k) = (1 - z^f +i-2)(^2^i+i-2)fc^ ^■<^2 

Fr{xi,n2+s = k) = {1 - aszf-^){aszf-^)'' is = l,...,p). 

(p) (p) 

Introduce the notations h- := fi- + n2 + p — i and 

In this setting, the joint probability that the matrix [xij]nix{n2+p) such that the sub-matrices 
[^i,j]nix{n2+s) {s = 0, ■ ■ ■ ,p) Correspond, under RSK, to pairs of tableaux of shape fx^'^^ has the 
explicit form 

i=l ^ ' 'h\'''' i<j i<j 

with 

^^-(n2+p)Ej=i 3 ^ 

U?ll\t;t)i U?ll\t;t)i 

m— 712— p— 1 ni n2 p ni 

X n n 11(1 - z^f^^-') n n(i - «^-*^"')- 

Z=l i=lj=l s=l i=l 

2.3 The Jacobi limit 

The Jacobi limit of tlie setting of Proposition[T]corresponds to each site of the ni x (n2+p) square 
lattice specifying a non-negative continuous exponential random variable with site dependent 
variance j < n2 given by 

Pr(xij e[y,y + dy]) = {i + j - 2 + 2a)e-^(^+^-2+2'^)dy, j < na 

FT{xi,n2+s€[y,y + dy]) = (i - 1 + a + a,)e-^(^-i+"+"^)d2;, {s = l,...,p). (2.13) 

This can be obtained from (|2.12|) by setting 

t = e'i/^, z = e-'^/^, as = e-''^/\ hf/L = xf (2.14) 

and taking the limit L oo. The joint probability ()2.12p . scaled by multiplying by L'^^+p^("2+p/2)^ 
has the following limiting form. 

Corollary 1. The PDF obtained by the limiting procedure {2.1^^ applied to \2.1S^) is equal to 
k^,^2^{{a^.a)e-'^^^^'^^^^'^f{e-'^^^^'^'-^^^^ > x^-^)) 

s=l 

X n(l-«"^'"r"" n n ie-^^'-e-'^') (2-15) 

i=l l<i<j<n2+p l<*<i<'^2 
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where 



K ( ^ .\ - nri^'"'' ^! rr + « + ^i) n r(2a + z + n2 - 1) 

"^•"^•^^^ ^"(nrir^/!)(nrir'/!),y nas+a) n r(2a+i-i) 

and 

> := > x^^) > • • • > > 4r+li > 4ts > 0). (2.16) 

Changing variables e~^J = yj'*'* the PDF (|2.15|) reads 

n2+p 712 P 

i=l j=l s=l 

ilj=l Vyj y l<i<j<n2+p \<i<j<n2 

where 

x(y(^) < := x(0 < yS^) < yt'^ < ■ ■ ■< ytli < vLts < l)- (2-18) 

In the special case Us = a — j {s = 1, . . . ,p) this simplifies to the functional form 

n2+p p 
i=l l<i<j<n2+p l<i<j<™2 s=l 

with C the normalization (C will be used generally below for this purpose and so its explicit 
value may vary from equation to equation) and w{y) = — y)^ for certain a and /3. The latter 
is the Jacobi weight, which in a functional form involving Vandermonde products is typical of 
a PDF arising in random matrix theory. Indeed, this functional form for each of the Gaussian, 
Laguerre and Jacobi weights can be obtained as eigenvalue PDFs. 

Before turning to such random matrix interpretations, it should be pointed out that in the 
setting of Proposition [U choosing each site of the ni x (n2 + p) square lattice according to 
continuous exponential random variables 

PT{xi,je[y,y + dy]) = e-ydy (2.20) 

leads to ()2.19p with the particular Laguerre weight w{y) = y"-i~("2+p)g-s/_ Further, rescaling 

the variables yj^'* i-^ ni(l + y^^ ^J2jn\) therein, and taking ni ^ oo gives ()2.19p back with the 

Gaussian weight wi^y) = . In the case p = 1 both of these facts are explicitly demonstrated in 
\14:\ Props. 4&5]; the extension to general p involves straightforward limiting procedures applied 
to (|2171) . 

3 Joint eigenvalue PDF for some nested sequences of random 
matrices 

3.1 Gaussian unitary ensemble 

By definition, matrices Mjv from the N x N GUE satisfy the recurrence 



+1 



Mat w 
a 



(3.1) 
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where a ~ N[0, 1/^/2] and each component Wj of w has distribution Wj ~ N[0, 1/2] + iN[0, 1/2]. 
Further, with Un the unitary matrix which diagonahzes Mat, and thus Mat = UnD^u\^^ where 
D]\j \s the diagonal matrix of the eigenvalues of Mjy, one has 



(3.2) 



This bordered form is the key to studying the joint distribution of the eigenvalues of the sequence 
of GUE matrices {Mj}j=i^2,...! or equivalently that of the sequence of principal minors of a single 
infinite GUE matrix Pllil lll]. 

In particular, it follows from (j3.2p that the characteristic polynomials pn{X)^pn+i{X) for 
Mjv,MAr_|_i are related by 





■ 




" Mn w ' 




' Un 


■ 


t 


Dn w 


0^ 


1 




w"^ a 




0^ 


1 




w'^ a 



A 



N 

E 



A- A 



(TV) 



where {A^^^} denotes the eigenvalues of Mn assumed ordered 



< A 



(iV) 



With {A'^''} regarded as given the PDF for the zeros of this random rational function, and thus 
the PDF for the distribution of the eigenvalues {A^^"*^^^} of Mn+i, can be computed to be equal 
to 



Af+l/,(Af+l) 



niV+i/ \ u 



A 



(N+l) 



AN) .{N) 



where xi^^'^^"^"^ < A^^-*) is specified as in (j2.18p except that the first and last inequalities are 
replaced by — oo < and < oo respectively. From this result, together with the fact that the 
eigenvalue PDF for N x N GUE matrices is proportional to 

1=1 i<j<k<N 



x(A(^+i) < aW) 



it follows that the joint eigenvalue PDF for the sequence of GUE matrices {Mn}n=N,...,N+p is 
given by (|2.19|) with n2 = N and w{y) = e~^^, and the definition (j2.16p of xiu^^^ < y^'^~^^) 
modified appropriately. 



3.2 Laguerre unitary ensemble 

Matrices A(„) from the N x N LUE with parameter a = n — N are constructed from n x N 
rectangular complex Gaussian matrices with entries N[0, l/\/2] + iN[0, l/\/2] according to 
^(„) = Such matrices satisfy the recurrence 



A 



(n+l) 



^(n) + Xx\ 



A 



(0) 



[0] 



AfxAf 



(3.3) 



where x is an x 1 column vector of complex Gaussians. Our interest is in the case n < of 
for which there are N — n zero eigenvalues. Then, after making use too of the invariance of 
xx^ by a unitary similarity transformation, the recursion (|3.3p can be written in the equivalent 
form 

= diag(ai, . . . , a^, 0, . . . , 0) + ff^ (3.4) 

N-n 
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where {ai}j=i^,.,^n are the non-zero eigenvalues of assumed ordered so that < ai < • • • < 
On- From (|3.4p it follows that the corresponding characteristic polynomials are such that 

det(Aljv _ ^ " [xjf Ejln+ikjf , . 

det(AlAr - ^(„)) ~ p^A-Oj A ' ^ ' ' 

The conditional PDF for the zeros of this random rational function, and thus the conditional 
PDF for the eigenvalues of can be computed exactly as [T^ [TT] 

n+1 n ^ "rr"+i^'\ _ \ ^ 

i=i j=iOi lli<iia» fljJ 

Recalling that the non-zero eigenvalue PDF for the matrices is proportional to 

n n 
j=l i<k 

it follows that the joint eigenvalue PDF for the sequence of LUE matrices {A(^j.)}r=n....,n+p with 
n + p < N is given by ()2.19p with n2 = n, w{y) = y^~("'+p)e~^. 



3.3 Corank 1 random projections 

Let An be an n X n matrix with eigenvalues ai < a2 < ■ ■ ■ < an, and let x be an n x 1 random 
complex Gaussian normalized column vector. The matrix 

Mn := n„A„n„, Hn ■■= In " xx^ (3-6) 



then represents a corank 1 random projection of An- We know from [T3] that in general M„ has 
a single zero eigenvalue, while the non-zero eigenvalues {\j}j=i,...^n-i have the conditional PDF 



, ni<j'(A.-A,) 



(n - i), ^^<^ ;"^ '"^' i^a < A). (3.7) 



Introduce a nested sequence of matrices {Ai}i=i,...,n by setting equal to the diagonal 

matrix formed from the non-zero eigenvalues of M„, computing M„„i according to (j3.6p . setting 
An_2 equal to the diagonal matrix formed from the the non-zero eigenvalues of M„_i, and 
repeating. With the eigenvalues of An having PDF proportional to 

n 

Hwiai) J] {ak-ajf (3.8) 

1=1 l<j<k<n 

it follows immediately from (13. 7p that the joint PDF for {Ai}i=n-p,...,n is given by ()2.19p with 
n2 + p = n. 

We remark (see e.g. [^) that the Laguerre and Jacobi cases can be realized by the matrix 
structure An = XnXn for certain matrices (recall the discussion of the previous subsection for 
the Laguerre case). It follows that the above construction is then equivalent to applying a 
sequence of corank 1 projections directly to the matrix X„. 
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4 Correlation for the classical projection process 
4.1 Approach via a general formula 



In studying correlations associated with (j2.19p . the most general case occurs when n2 = 1. After 
then changing notation by setting p = N — 1, y'f^ ^ x^p^^li^ (|2.19p reads 

N N-1 

^n^(^S''^) n {^T-A''^)\{x{x^'^'^>x^'^ (4.i) 

1=1 l<j<k<N s=l 

Here x{^^^~^^^ > x^*)) is defined as 

(cf. (j2.16p ') and w{x) involves a factor Xo<a;<i- Hereafter we consider more general cases in 
which w{x) does not always involve such a factor (see (14. 7|) below). Of course the product of 
differences can be written in terms of the Vandermonde determinant, giving 

N-l 

n '^(^l^^) n ^^f^ ~ 4^^) °^ det[w;(xj.^^)p7v-j(4^^)]j,fc=i,...,iv, 

1=1 i<j<k<N 

with {pj(2;)}j=o,...,Ar-i a set of arbitrary polynomials, Pj{x) of degree j. Furthermore, we know 
from [13, Lemma 1] that 

^(^(.+1) > ^(«)) = det[x{xf^'^ > xi'^)],^k=i,...,s+i (4.2) 
where 2;^^^ := —00. Consequently ()4.ip can be written in the form 
1 

- n det[,^(x;.^), 4^+^))],,, =i„,.,,+idet[*j^_,.(xf ))],,, =i„„,^ (4.3) 



s=l 



with 



Hx,y) ■=Xy>x, ^fix):=w{x)pj{x). (4.4) 

The general structure (|4.3p is precisely that for which the correlations have been determined 
in the recent work [U Lemma 3.4]. To apply this result, with {a*b){x, y) := a{x, z)b{z, y) dz, 
it is necessary to compute the quantities 



'"2) (x,y) := (</.*••• *</>)(x,y), m < 



"2 



n2— ni times 

(for ni > 722, (/)("i'"2)(x, y) := 0), and 

C-,(x) := (0("'^) * *]^-,)(^) (1 < n < iV, J = 1, . . . , iV). 
Use of shows 



(n2-ni-l)! 
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(with the convention that l/{—p)l = for p £ Z^, this vanishes for ni > ^2) and 



To proceed further, we choose w{y) to be one of the classical weight functions 

, Gaussian 
w{y) = { y'^e-yxy>o, Laguerre 
y"-{l - yfxo<y<i, Jacobi. 



(4.6) 



(4.7) 



We further choose Pj{y) to be proportional to the corresponding orthogonal polynomials, as 
specified by their Rodrigues formulas 



Vj{y) 



1 

ejw{y) dyi 



w{y){Q{y)y 



Hj{y), Gaussian 
L^j^\y)^ Laguerre 
_ Pj"'*^(l-2y), Jacobi 



(4.8) 



with the quantities ej and Q{y) defined in the various cases by the pairs 

((— I)-',!), Gaussian 
(ei,Q(y)) = < (j!,y), Laguerre 
^ (2^j!,y(l-y)), Jacobi. 

Substituting (j4.8p in (j4.6p and integrating by parts shows that for j > (n 7^ A^) 

w{x)Hj{x), 



(4.9) 



e-N-n+j 



w{x)\a^a + N-n.P) 



(a+N-n 



b+N-n) 



Gaussian 
Laguerre (4-10) 



(1 — 2x), Jacobi, 



while for j < 



_l'^N-n+j I 
eN-n+j {-j-iy. 



(y - x)-^-it«(y)(Q(y))^-"+^' dy. (4.11) 



As further required by ^ Lemma 3.4], one introduces the polynomials {<l'"(x)}j=o,...,n-i5 n 
1 , . . . , — 1 by the orthogonality requirement 



/oo 
^]{x)n{x) dx = 6j,k. 
-00 



From (|4.1U|) we see that 



$"(x) = (-1) 



where 



1 p 

ATj I ai-^a + JV-n J 



r(i+a+i) 



{a+N-n,b+N-n) 



Gaussian 
Laguerre 
(1 — 2x), Jacobi, 



(4.12) 



Gaussian 
Laguerre 



r(j+i) ' 
r(j+a+i)r(i+&+i) , , . 

j\{2j+a+b+l)r{j+a+b+l) ' Jd'-UUl 



(4.13) 
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is the normalization associated with each polynomial respectively. 

A crucial feature exhibited by (j4.12p is that $q(x) is a constant. Now, Assumption (A) of 



[H Lemma 3.4] requires that ^q{x) oc (j){x^^-^, x). Recalling from below (j4.2p that 



-OO, 



we see from the first definition in (j4.4p that indeed (j){x^_^-^,x) is similarly a constant. With 
Assumption (A) satisfied, [H Eq. (3.25)] gives that the correlation between eigenvalues of species 
Sj at positions yj {j = 1, . . . ,r) has the determinant form 



p(.{(.Sj,yj)} 



det[K{sj,yj;sk,yk)]j,k= 



:l,...,r ) 



(4.14) 



with the kernel K given in terms of the quantities (j)^'^'^''^^\x,y), ^'"(x), $"(x) specified above 
according to 



Sk 



K{sj,yj;sk,yk) 



{yj,yk) + ^KU{yj)'^:U{yk). 



(4.15) 



1=1 



These findings can be summarized in the following statement. 



Proposition 2. Consider the joint PDF 1^. j| ), with w{y) one of the three classical weights ^J^- 
Specify 5y g^. e.^Q(y) g^. ^n(^) g^. ^n(^) 

A/^- 62/ (^TT^p. /n terms of these quantities, the general r -point correlation is given by ^-1^^ with 
kernel ( [^.i5[ j. 



4.2 Direct approach 

Here the method of [26j will be used to reclaim Proposition [2j The starting point is (|4.3p . 
rewritten in the form 



1 

C 



N 



JJdet 



s=2 



1 



(Af-s)x(Af-s) 







sx(Af-s) 



^(N~s)xs 
{x^^~'^\x^u'') - Ksix*"^'^^ 



)] J = l,...,s-1 
fc = l,...,s 



(4.16) 



[l]fc=l,...,s 

X det[V'7_i(4''0]i,fc=i,...,Ar 

so that all determinants are of the same dimension. The auxiliary function Ks{x) is arbitrary, 
the value of the determinant being independent of Ks{x). 

In the notation for w{y)^ Pj{y) Mj as defined by (|4.7p . (j4.8p and (|4.13p . introduce the 
superscript (N—n) to indicate that a 1— > a+n (Laguerre case), a a+n, b ^ b+n (Jacobi case). 
In the Gaussian case the superscript has no effect. With this meaning understood, expanding 



{x,y) = Xy>x in terms of {pj''^(y)}j=o,i,... gives 



fc=0 -^k 



(s) 



' w^'\t)pi'\t)dt. 



(4.17) 



Separating off the A; = term, making use of the Rodrigues formula (j4.8p and integrating by 
parts reduces this to 



4>{x,y) 



u;W(t)dt-^(-^)(x)^ p[^-'\x)pi%,{y). (4.18) 

K -^^^ fc=o ek+l^fUl 



Recalling from (j4.16p that k,s{x) is to be subtracted from (j){x,y), the formula ()4.18p suggests 
choosing 

1 



Ks{x) 



Ml 



W^'\t)dt. 



(4.19) 



11 



Making this choice, and with the notation 



(4.20) 



(note that {^l^\x)}k=o,i,... is a set of orthonormal functions) we obtain for (j4.16p the expression 



N 



c 



JJdet 



s=2 



1 



(7V-s)x(7V-s) 







sx(N-s) 



^(N-s)xs 
[l]fc=l,...,s 



= 1,...,3-1 
fc = l,...,S 



xdet[7?5.!^J(4^^)]j-fc=i, 



where, with 

the quantity (ps is specified by 



(4.21) 
(4.22) 



oo (s— 1) 

fc=o 7fc+i 

and for convenience the final row in the bottom right block of the first determinant has been 
moved to the first row. 

Our next step is to introduce 



Sj^i-i, otherwise 



and in terms of this to define 



N-l {s) 

A(s,t) _ Sr^ ^k+s-N (s) (t) 
j,l '~ (t) ''^k+s~N,j+s-N''lk+t~N,l+t~N 

k=0 Ik+t-N 



CO (s) 

E ^k+s-N (s) (t) 
(t) ''lk+s-N,j+s-N''lk+t-N,l+t-N- 
k=0 lk+t~N 



We can write (|4.2ip in terms of {Ajl }, {G - / } so that it reads 

TV 

M)l 



1 

C 



s=2 



c 



det 



^(iV,l) ^(7V,2) ^(iV,3) ^(iV,4) 

-G(1'^) 

-G(2.3) _G(2.4) 





j^{N,N) 
_QiN-l,N) 



(4.23) 



12 



where A^^'*) := [Ajf']j^i=i^,„^N, G^'''*) = [Gjf']j,i=i,,„^N- With a^^'*) the iV x iV matrix such that 
Q,('5:t)^(*>") = j{is,u)^ multiply row 1 by a^^~^'^^ and add to row j {j = 2, . . . ,N) to rewrite this 
as 

j^{N,l) j^{N,2) ^(7V,3) ^(7V,4) ... j^(N,N) 



c 



det 



^(1,1) ^(1,2) ^(1,3) ^(1,4) 

^(2,1) ^(2,2) ^(2,3) ^(2,4) 

^(3,1) ^(3,2) ^(3,3) ^(3,4) 

^(7V-1,1) ^(JV-1,2) ^(iV-1,3) ^(JV-1,4) 



5{3,iV) 
^(7V-1,JV) 



where := ^(*'*) — G^^'^\ Moving the first block-row to the final block-row gives the 

structured formula 



idet[F(' 



s,t) 



(4.24) 



Moreover, from the definition of A^*'*) and G^*'*-* we observe that for s > t 



j<N-soi:l<N-s 



while for s < t 



0, 



j<N-soTl<N-t. 



This allows the dimension of the block matrix in (j4.24p to be reduced, giving for the joint PDF 
(SH), P say, 



C 



det[/ 



{s,ty 



s,t=l,...,N 



where x^-'^ = {x^ , . . . , Xj ) and is the s x t matrix with entries 



j-s+N,l-t+N 



2^k=i jt) Vs-k\^j )\~k\^i A 



s > t 



E 



fc=-oo (t) Vs-k\^j l^t-k^^l h S <.l. 

't-k 



From the orthonormality of {i]k\^)} follows from (j4.25p that 



As,t) At,u) J it) 

Jj,l Jl,m "'■^l 



S>t>UOIS<t<U 



As,u) 
J j,m ' 

0, otherwise. 



(4.25) 



(4.26) 



(4.27) 



We seek to use the form (j4.25p . together with the property ()4.27p . to compute the correlation 
between eigenvalues of species sj at positions yj (j = l,...,r). For this purpose, we group 
together the eigenvalues of distinct species in the correlation. Thus if the distinct species are 
si, . . . , Sf, with §1, . . . ,Sf € {si, . . . , Sr}, we write the positions being observed in species s as 



{x^^\ . . . ,Xn^) (1 < rig < s). The correlation relating to is specified in 



terms of the PDF p by 



P{{x^'^%=l,...,r) 



xp{x'^^\...,x(''>). 



N 

n 



s = l 
s^{si,....Sf} 



dx 



(4.28) 
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Because of the structure ()4.25p and the orthogonahty relation (j4.27p . these integrals can all 
be computed by performing a Laplace expansion of the determinant (see e.g. [MIIH]) to give 



L k — l,...,ni^ 



a,f3=l,...,f 



In the notation of (j4.14p this can equivalently be written 



p{{{sj,x,)}j=i_r) = det[f^^^''\k=l,...,r (4.29) 

(s ■) 

(note that the superscript on Xj is now redundant, and hence has been omitted) so to complete 
our task of rederiving (j4.14p it is sufficient to show that 

^'^'' = SfSy^^^^'"^'^^'"^^ (4-30) 

(the corresponding determinant is independent of the function a{s,x)). 

To verify (j4.30p we begin by recalling (j4.20p and (j4.22p to see from the definition (j4.26p that 
for s >t 



fj,i = 'i^jW 'i^i)) 2^—- -Jt) (4-31) 

fc=i ""'-^ ^t-k 



and for s <t 

fc=-oo * ^t-k 

On the other hand, it follows from (liTSll and (14.101) that for j > (n / iV) 



C(x) = (-1)^-"— ^u;(")(x)p;")(x), 



while according to (|4.12p 



^]{x) = (_i)-n£^Xz!i±l^. (4.33) 
Recalling too that (j)^^'^\x, y) = for s >t we then see from the definition ()4.15p that for s > t 



Kis,xy,t,x,) = (-l)-V-)(x,);^^ "-^^-^^;"*-^^-^^ (4.34) 



k=i ^'-^ K^-k 



r(t) 

k=l ^'^ 

For s <t 



<j)^''^\xj,Xl) = J^^J—^Xx,>x,{xi - XjY ' \ 

Analogous to (j4.17p . we can expand in terms of {p^*^(y)}, 

k=0 4' 



^ i,_n J^i. 
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Proceeding now as in the derivation of ()4.18p gives 

k=-oo ^t-k^^t-k 

+ E a !"T iv''[% r-(*-)(.)(.-x)— dn. (4.35) 
But from (|4.1ip . for j < 

m{x) = ^ — '- / {u - x)-^-^w^''~^Hu) du 

CN-n+j {-j -ly- Jx 

so after making use too of ()4.33p we have 

t 

p=s+l 

= E' - "y. (4.36) 

k=0 ^ f^kJ^k 

Thus adding (j4.36p to minus (j4.35p cancels the final line in the latter. It remains to add to 
minus (|4.35p the quantity 

EIS I / \ I T\S-t Sf \Sr^ ^S-kP^S-ki^)p\-kiy) orrN 

P=i fc=i "^^-^ ^t-p 

This cancels the corresponding terms in the first sum of minus (I4.35p . giving the result that for 
s <t 

K{s,x,;t,xi) = {-lY'^-'w^^Hx) ^ pi%{x)pfl,{y). (4.38) 

fe=-oo ^t-k'^t-k 

Comparing (fOT]) with (Oip . and (f02]l with (f08]l . we see that for general s,t 
thus verifying (I4.30p . 

5 Scaling limits 

It is well known (see e.g. [9", "8]) that the eigenvalue distributions for the joint PDF (j3.8p with the 
classical weights (14. 7p permit three distinct scalings as n ^ oo. These correspond to eigenvalues 
in the bulk of the spectrum, or in the neighbourhood of the spectrum edge. There are two 
distinct cases of the latter — the soft edge and the hard edge. The hard edge is characterized 
by the eigenvalue density being strictly zero on one side. This occurs for j; < in the Laguerre 
ensemble, and for both a; < and x > 1 in the Jacobi ensemble. In contrast, the neighbourhood 
of the largest eigenvalue of the Laguerre and Gaussian ensembles is such that the eigenvalue 
density is to leading order in n zero, but at higher order it is non-zero. This is referred to as a 
soft edge. 
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For the projection process (j4.ip with classical weights ()4.7p we again expect these same three 
distinct scalings, provided the difference between ranks of the matrices (or equivalently between 
labels of the species) is fixed. We find at the soft edge the correlations can be interpreted 
as though the eigenvalues of the diff'erent species coincide with the eigenvalues of species (A^). 
Thus they are fully determined by the Airy kernel (see (|5.3p below). That the eigenvalues of the 
different species should coincide at the soft edge is not at all surprising upon consideration of the 
joint PDF ()2.19p . Thus one observes that the lowest indexed species repel via a Vandermonde 
factor, with no restoring potential apart from the ordering constraint. Thus they will tend to 
cluster at the boundaries, which at the soft edge corresponds to the positions of the species (N). 

The situation in the bulk and at the hard edge is more delicate in that the correlations 
depend on the difference between labels of the species, even though this difference does not scale 
with N. However the correlations within a given species can be anticipated. With w{x) = x°'e~^ 
in (14. ip the marginal distribution of species — c is precisely the Laguerre unitary ensemble with 
a I— > a + c. Thus the hard edge correlations within species N — c must be the usual Bessel kernel 
correlations with a ^ a + c, which is indeed what we find. Because the bulk correlations for 
the Laguerre and Gaussian unitary ensembles are given by the sine kernel, all bulk correlations 
within a species will be specified by the usual sine kernel. As to be discussed below, for the 
bulk scaling it turns out that the full set of correlations are those known from the so called bead 
process [7J. 

Soft edge scaling, together with the difference between labels of the species scaling with N is 
of particular interest for its relevance to the queueing process of Baryshnikov [Ij, or equivalently 
a lattice version of the last passage percolation model of Hammersley (see e.g. |10j). To see this, 
let us revise some aspects of the theory relating to the latter topic. 

Thus, with each site (i,j) in the quadrant x Z"*" an exponential random variable xij of 
density 2e~^*, t > 0. Define the stochastic variable 

Z(m, n)=max Xij, (5-1) 

(1,1) u/r (r?i,n) 

where the sum is over all lattice paths in Z+ x Z+ which start at (1,1) and finish at {m,n) 
going either one lattice site up (u), or one lattice site to the right (r). It is well known that with 
Xn '■= l{n,n) the variables {xi, . . . ,Xn} have a joint PDF of the form (j3.8p with w{x) = 
\18\ . Hence, as n — > oo the corresponding distributions permit a soft edge scaling describing the 
scaled distribution of l{n,n). This scaling is fully described by the correlation function 

pl^^ivi, . . . ,y„) = det[K'°'\yj,yk)U=i,...,n (5.2) 

where 

.soft. s . Ai(x)Ai'(y) - Ai(y)Ai'(a;) 



i^-"(x,y) := 



X -y 



oo 



Ai(x + u)Ai(y + u) du (5.3) 







is the so called Airy kernel. It is also known [19j that the sequence of stochastic variables 
{l{n + k,n — A;)}fc=o,i,... permit a scaling to a state specified by the dynamical extension of the 
correlation (|5.2p . 

/3^°f((ri,yi),... ,(r„,yO) = det[K'"'\{Tj,yj),{Tk,yk))]j,k=i,...,n (5.4) 
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where 



/•oo 

A^^\x,y) := / e"™Ai(x + n)Ai(y + u) 
Jo 

A^^\x,y):=- e^™Ai(x + n)Ai(y + u) (in. (5.5) 

J —CXI 

This is the so cahed Airy process A2, which underhes the distribution of the largest eigenvalue 
in the scaled limit of the Dyson Brownian motion model of the GUE |25t . 

The significance of these facts with respect to the present study is that in ()2.19p with n2 = 0, 
rii = p = n, w{y) = e~^, we know from the sentence including (j2.20p that the variables {yj" 
(i.e. the largest eigenvalue for each species) coincide with the stochastic variables {l{n,n — 
A;)}jt=o,i....- By analogy with the behaviour of the stochastic variables {l{n + k,n — fc)}fc=o,i,..., 
one may anticipate that their distribution is controlled by the Airy process A2 ■ We will find that 
with the differences between the ranks (species) scaled to be of order n?^^ that this is indeed 
the case, and that the same effect holds for the soft edge in the Gaussian case. 

5.1 Fixed differences between species 
Soft edge and bulk scaling 

Explicit details will be worked out only in the Gaussian this case is typical; in particular 

the scaled correlations do not depend on the particular case they originated from as is typical 
in random matrix theory (a form of universality). In the N x N GUE the soft edge scaling is 
obtained by the change of variables 

Xi = Vm + ^ \ . (5.6) 
^/2iVi/6 ^ ' 

This has the effect of moving the origin to the neighbourhood of the largest eigenvalue, and 
scaling the distances so the inter-eigenvalue spacings in this neighbourhood are of order unity. 
The bulk scaling is obtained by the change of variables 

(5.7) 



/2iV' 



which makes the mean particle density in the neighbourhood of the origin unity. With the 
species differing from by a constant, 

Si = N-Ci, (5.8) 
we seek the limiting forms of the correlation (j4.14p for both the soft edge and bulk scalings. 
Proposition 3. For the soft edge scaling 



and with Si specified in terms of Ci by 115. 8\) . 



K{sj,yj; si,yi) ^ K'°'\Y„ Yi), (5.10) 



V2iVV6"^'^^'-^'-'-^^-.oo aNici,Yi) 
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where 7^^°^' is given by i5.3\) and aN{c,Y) := e ^^^''^^{2N) "^Z^. Consequently 

i™oo (7I^)'^(^("^'y^)^^=i'-'^) = det[K-f*(y„n)],,fc=i„„,,. (5.11) 
For the bulk scaling 

Vi = 4=, (5.12) 
and with Si specified in terms of Ci by 15. 8\) . 

^i^(.,,,,;.,,0,^^^i?((c„y.),(Q,y.)), (5.13) 

where bN{c) := {2N)~'^/'^ and 

B{{t^,x), {Ty,y)) := < 



/ s^y cos(7rs(j; - y) + 7r(T^ - Ty)/2) ds, Ty > 
Jo 

gTy-T^ cos(7rs(x - y) + 'k{tx - Ty)/2) ds, Ty < t^^ 



(5.14) 



Consequently 

J™ (^)'/o({(5i,yi)}j=i,...,r) = det[Bi{c,,Y,), (Q,yO)ki=i,...,r. (5.15) 

Proof. Substituting the Gaussian case of (j4.7P " (j4.9p and (|4.13p in (j4.34p shows that for Sj > si 
{Cj <ci) 

K{sj,yf,si,yi) = 2^1-^(1^ _ j^y H,^^k{yj)Hs,-k{yi)- (5.16) 

Our strategy is to use appropriate expansions of the Hermite polynomials, corresponding to the 
different scalings, to simplify the summation. 

Consider first the soft edge scaling. With x related to X by (|5.6p we have the uniform large 
N expansion ^30j 



e--'/2i/;v(x)=vrV42^/2+V4(iv!)i/2^r-i/i2|^Ai(X) + 0(iV-2/3)| [^|^^'')' ^Jq). (5.17) 
We rewrite this to read 



and then substitute in (j5.16p to obtain 

K{s„yy,sum) ~ e-^^^^«-^^)2-fe-')/22V2,-V6 

X E ( I^Ic^Ifc)' ) '^'^'^^^ + fc/ivV3)Ai(y, + (5.19) 

The leading order contribution to the sum in (|5.19p comes from terms o{k^^^). Noting that 
then (A^ — Cj — k)\/{N — ci — k)] ~ N^^~'^^ the sum can be recognised as the Riemann sum 
approximation to the integral form of K'^°^^ in ()5.3p . and ()5.10p in the case Cj < q follows. 
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We turn our attention next to analyzing (j5.16p in the bulk scaling limit. For this we use the 
uniform asymptotic expansion 

^M^±^e-^'I^Hn{x) = cos(^/2^x - nvr/2) + 0{n-^'^), 
r(n + 1) 

and a simple trigonometric identity to deduce that 



1 

K{sj,yj;si,yi) ~ — = ^ 



1 {N - Cj - k)\ 

2\/^ ^ ({N - c,- - k)/2)\{{N -ci- k)/2)\ 



X cos (vry ^^j^iYj - Yi) + - ci] 

Here the main contribution to the sum comes from (N — k)/N = 0(1). Expanding the ratio of 
factorials in this setting gives 



N 



/ N — k TT \ 

X COS (vr^ ^^(y, - Yi) + -{cj - q)) . (5.20) 

Recognizing the sum as the Riemann sum approximation to a definite integral in the variable 
[N — k)/N = t, then changing variables t = in the definite integral gives (j5.13p for Cj < c/. 

It remains to study the case Sj < si {cj > ci) for both the hard and soft scalings. For the 
soft edge scaling it turns out that the form (|4.38p is not appropriate. Instead we make use of 
the form ()4.15p . which recalling ()4.36p and (I4.37P can be written 

K{s,x;t,y) = -4>'~'''\x,y)+ a - k - IV^' W / " ^)''^''~'^^'"'^ W 

fc=0 ^ ^ '' ^k^k 

+(-l)-V^)(x) y 2fzt P's-ki^^pf-k(y) _ (5 21) 



Considering this as the sum of three terms, the first two do not contribute in the scaling 
(15.81). and so 



' e,^_kP''s--kiyj)P^-kiyi) 



K{s„yr,si,m) ~ (-l)^.-^u;(^.)(y,.)E7^- 

This is precisely the expression ()4.34p . except that the upper terminal is sj instead of s/. Recall- 
ing the working below (|5.16p . this detail does not affect the leading asymptotic form, so (|5.10p 
applies for both Cj < Ck and cj > Ck- 

In distinction to the strategy required at the soft edge for Sj < si {cj > q) , to analyze the 
bulk scaling in this case the form (j4.38p is well suited. As the only difference between (j4.38p 
and (14.341) is in the range of summation, the working leading to (j5.20p again applies, so this 
asymptotic formula remains valid but with k summed from —oo to 0. Crucially, because Cj > ci 
this sum is convergent (it is in relation to this requirement that an analogous approach to the 
soft edge scaling breaks down), and is furthermore a Riemann sum approximation to the same 
definite integral as found for the case Cj < q, but on (— oo,0] instead of [0, 1], hence implying 
the second formula in ()5.13p . □ 
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5.2 Hard edge scaling 

Hard edge scaling is possible for both the Laguerre and Jacobi cases; here the details will be 
given in the Laguerre case only, as the limiting correlations are the same in both cases. For the 
N X N LUE the hard edge scaling results from the change of variables 

X. = §, (5.22) 

which makes the inter-eigenvalue spacings in the neighbourhood of the hard edge a; = of order 
unity. We seek the limiting correlations with the scaling (j5.22p and the species specified by ([57 



Proposition 4. For the hard edge scaling 15.22\) and with Si specified in terms of Ci by 15. 8\) 

■^K{sj,Xj;si,xi) ~ ^^^^^y^H{cj,Xj;ci,Xi), (5.23) 



where hN{c) := (2A^) ^ and 

' 1 

i?((rx.,x), (ry,y)) := < 



Consequently 



S^^y-^^'^l^Ja+rA{sxY^^)Ja+r,{{syY'^)ds, Ty > 

. ^ .oo (5.24) 

/ s^^y-^-^/\ja+rA{sxf'^)Ja+r,{{syf'^)ds, Ty < T,. 



J™ (tK?) P{r){{{Sj,Xj)]j=l,...,r) = '^ei[H{cj,Xf,CuXl)]j^l=i^,,,^N- (5.25) 
N — ^oo V 4iV / 



Proof. Substituting the explicit form of the Laguerre case of the quantities in (j4.34p gives 



N-ci 



KiN-c„xr,N-c,,xO=x;^^^e-^^ ^ _l ^ ^^"^^ L^,(x,)L^,{xO, (5.26) 

valid for q > cj. As scaled according to ()5.22p , it is appropriate to make use of the 

uniform asymptotic expansion [32] 



-x/2 al2Tai \ al2i i^i \\I2\ , f X^/^0(n"/2 ^Z"^), CU ^ < X < UJ 

e ^/^x'^^^L'^ix) =n''/^Ja{2{nxy'^) + ■{ a/2+2r^fa\ ^ -1 

1^ x"'^+^0(n"), < X < cn \ 

Using this, and expanding the ratio of gamma functions with (N — k)/N = 0(1), we deduce 

K{N - Cj,Xj]N - ci,xi) 

, .r , ^V^^-^A^^'""'^^^ ffN-k \1/2n ffN-k \l/2 



k=l 



(5.27) 



This is a Riemann sum, and the result (j5.23p in the case q > cj follows. 

The expression (I5.26P is also valid for ci < cj, provided the summation in now made over 
k £ Z<o. Following the above working through again gives (j5.27p . but with the summation over 
k G Z<o. Because ci < Cj the sum is convergent, and its leading form given by the definite 
integral made explicit in (I5.24p . □ 
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5.3 Soft edge scaling with difference between species 0{N'^^'^) 

As anticipated from the viewpoint of last passage percolation, the soft edge scaling permits well 
defined correlations with the species separated by 0{N'^^^). The details can be worked out for 
both the Gaussian and Laguerre cases, although the limiting correlations correspond to the Airy 
process A2 and so are independent of the particular case. 

Proposition 5. In the Gaussian case, scale Si according to 

Si = N + IciN"^/^ , (5.28) 

and scale yi according to 



y, = (25,)'/' + -^. (5.29) 



For large N 



where are given by ^ and aN{c,Y) := e-^'^'^(2iV)^^'^'e^'^''='e-2='/3. Conse- 

quently 



lim ^ 



N-ooi^/^mJ^) P({(«.'yj)}j=i,...,'-) = det[i^^°^*((-c,,X,),(-Cfc,Xfc))],-fc=i,...,,. (5.31) 

In the Laguerre case, scale Si according to 

Si = N-Si, Si := 2ci(2iV)2/3 (5.32) 

and scale y\^''^ according to 

y\''^ = Asi + 2{a + N- Si) + 2(2iV)V3y.. (5.33) 

For large N 

2{2N) I K{sj,yj;si,yi) ~ ^TJT^ \ A) \^ ^ (5-34) 

Pn{ci,Yi} [ Al/_^^{Yj,Yi), ci<Cj 

with Pn{c,Y) = e-(2^)'/='yAr-c(2iV)2/3g2(27V)V3,2^8c3/3, Consequently 

lim (2{2Nf/^)y p{{{sj,y,)},=i_r) = det[K''''\{c,,Yj),{ck^^^^ (5.35) 

Proof. The derivation is very similar in both cases, so we'll be content with presenting the 
details in the Laguerre case only. Reading off from (j5.26|) we have 

KiN - s,,y,- N - sum) = e'^^ ^ -t + a + ^ ^S^-S-^^^^^^S^^Ul^O- (5-36) 

fc=i 

This formula is valid for Sj < s/ (for sj > s/ the RHS is to be modified by multiplying by —1 and 
changing the summation terminals to A; S Zf<o; this modification does not change the working 
below in any essential way, and so will not be considered explicitly). 
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We seek the asymptotic form of (|5.26p upon the scahngs (|5.33p and (j5.32p . Adapting the 
strategy of the proof of Proposition [3l our chief tool is the uniform asymptotic expansion [23] 

^a/2g-x/2^a(^) = {-1^ (271)-^'"^ ^/ {n + a)\ / n\ 

Ai,x, + o,„-/3){o<-;)' m ) i,,r, 

where 

x = An + 2a + 2{2nf/^X (5.38) 

(this form allows for a = o(n); the classical Plancherel-Rotach type formula given in e.g. [32] 
requires a to be fixed and correspondingly has ^ {n + a)\/n\ replaced by n"/^). Use of this 
formula, rewritten to read 

^a/2g-x/2^a_^(^) = (- 1)""'= (2n) V(™ -k + a)!/(n - k)\ 



with n = N — Si shows that for large N 



/(N—~ —kV\^^^ 
x(2iV)-2/3 ^ j Ai(y, + 2fc/(2iV)V3)Ai(y, + 2fc/(2iV)i 



(cf. (|5.19p ). The leading order contribution to the summation comes from k of order N^^^. Using 
this fact, noting from Stirling's formula that for large s 

{s-kyj 

and using this formula with s = N , ki = k + Si {i = j, I) the sum is recognised as the Riemann 
sum approximation to A^^^ as defined in (|5.5p . implying the result ()5.34p in the case q > Cj. □ 

6 Discussion 

As pointed out to us by A. Borodin, replacing (|5.14p by 



r 1 



1 



^ -^-1 (6.1) 



1 



-1,1] 



{isYy-^-e''^^''-y^ ds, 



leaves the determinant (|5.15p unchanged. Further changing scale x i— > x/TT,y i— > y/vr removes 
the vr's from the exponents, and multiplies each integrand by a factor of l/vr to account for 
the corresponding scaling of the correlation function. The significance of this form is that it is 
identical to the 7 = case of the correlation kernel for the so-called bead model pf, Thm. 2]. 
In fact the bead model was already known to be closely related to the GUE minor process [71 
Section 4.1]. The form ()6.ip can also be obtained as a limit of the incomplete beta kernel of 
Okounkov and Reshetikhin [28^ Section 3.1.7] (write the parameter 2; as z = 1 + ia, change 
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variables w = 1 + ias, rescale the space variable / by a^^ and take a — > 0). The recent work [16] 
obtains the incomplete beta kernel in the context of a study of random lozenge-tilings. Further 
the kernels of O Thm. 4.4] permit degeneracies to (j6.ip . 

Another discussion point is in relation to consistency between the present results, and results 
from ^13j. In [13j the correlations for the p.d.f. 

1 " 

_^e~^^.+y.)/2^Ai,.,-y,)/2 -Q ^^^^^^ (6.2) 

were computed, along with the scaled limits at the soft and hard edges, and in the bulk. The 
p.d.f. (j6.2p with A = —1, is identical to the p.d.f. (j2.19p with w{y) = e~^, n2 = n, p = 1 and 
Un+i — 0- Setting yl^Hi = would not be expected to alter the soft edge and bulk scaling limits, 
so it should be that the scaled correlations in [13] contain as special cases the results ()5.1ip and 
(|5T3]) for \cj -q| = 0,1. 

To see that this is indeed the case, we recall from [13j that with 

J V^a/vr, bulk 
^1 a/2(2n)i/3, soft edge, ^^'"^^ 

the scaled correlation for ki variables of species type x, and k2 variables of species type y was 
calculated to equal 

P(fci,fc2)(^i' • • • ^^ki;Yi, . . . jYfcj) 



det 



(6.4) 



where 



^scaled (y^y/) = K^calcd^y^y/) 

^scaled^y^ X) = -e"(^-^)xX>Y + e"^/^ / g-ai,/2^scaled(^^ y) dv 



Kf'^'iX, X') = _e"(^-^')/2^ j g-a./2^scaled(^^ ^) _ (g 5) 

In the soft edge case K^^^^"'^ = K^oft specified by ([53]), while in the bulk K'"'''^'"^ = K^""^^ 
where 

.bulk.x. _ sin7r(X -Y) _ 



^buik(^^ y) = ; / = / cos 7r(x - y)t (6.6) 

We see from (16. 3p that A = — 1 corresponds to a = in the bulk, and a — oo at the soft 
edge. We see from ()6.5p that for a — oo 

K'^f\Y,X) -K'"^\Y,X), Klf{X,Y) -K'°^\X,Y), 

a 2 



Kl°J\X,X') ~ K'''^\X,X'). 
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When substituted in (j6.4p the factors —2/a, —a/2 cancel, and so agreement with (jS.lip is found. 
Further, setting a = in (|6.5p and recahing (j6.6p gives 



a=0 



a=0 



-XX>Y + 



X 



sin7r(v — Y) 
tt{v - Y) 



dv, 



TT I tsimriX -Y)tdt, ^^^^^"'^(X, X') 



a=0 



sin7r(X -X') 
7r(X-X') ■ 



A simple calculation shows that the first of these can be rewritten 



k^o'Hy,x) 



a=0 



sinvrf (X — Y) 



dv. 



TTV 



With this we obtain agreement with (j5.15p in the case \cj — q| = 0, 1, as expected. 
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Appendix 

Since the completion of this work, Borodin and Peche have posted a work [6] on the arXiv which, 
amongst other results, establishes our Proposition O The strategy used is, at a technical level, 
quite different to that adopted here. 

In this appendix we concern ourselves with another aspect of the work [6j, relating to a 
generalization of our ()3.3p . 

^(n+l) = + ^(n)4n)' ^(0) = Mpxp (^-1) 

where is a p x 1 column vector of complex Gaussians with entries such that the modulus 
of the i-th component has distribution r[l, l/('7rj + 7r„)]. As in Section [3.21 the point of interest 
is in the joint eigenvalue PDF for {^i, . . . , Ap}. This is not computed directly, but as in the 
discussion around p.20p , it is noted that the directed percolation in the pxp square which each 
lattice site containing an exponential random variable of density (vTj + Ttj)e~^^^~^'"i^^ has 
the distribution of the stochastic variable l{p,p) equal to that of the distribution of the largest 
eigenvalue of ^(p). For the percolation problem, the joint distribution of {/(j, p)}j=i,...,p can 

be calculated, leading to a joint PDF for the p species of variables {x^j^^}, {s = 1, . . . ,p) with 
j = 1, . . . , s proportional to 

det[e--4'''],,=i,...,/ndet[e-^'=-(^^''-^'''\,(.+i)>,,J.. ^ ^ 6"-^^^ (A.2) 



fc = l 



The question is asked as to whether this joint PDF is in fact the joint PDF for the eigenvalues 
of the matrices ^(s), s = 1, . . . ,p. 

In fact the working from |14^ Section 5] allows this question to be answered in the affirmative 
in the limit vrj — > c, (i = 1, . . . Thus it follows from \14\ Corollary 3] that the condition PDF 
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for the eigenvalues {aj}j=i^,,,^n of is proportional to 

n+1 n ^ Tr"+1/ 



n ^r^'^"''^ n ^e-(^+-)(^"i^^-S;-^^)^4^^x(A < a). (A.3) 



Let us now write Xj i— > A^"^^\ aj i— > Aj"^ The sort joint PDF is the product from n = 
1, . . . ,p — 1 of the conditional PDFs (jA.Sh . multiplied by the PDF in the case n = 1, which is 
proportional to (^x^^'^)P^^e^^'^^'^''-^^^'^ ^ . This gives (]A.2p with the first determinant therein replaced 
by nf=i e~'^^' YVi<j(.^^j''' ~ ^i^^)' thus verifying ()A.2p in the case that —>■ c, {i = I, . . . ,p). 
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